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Beta-lactams are used as major therapeutic agents against a number of infectious agents. Due to widespread 
use of p-lactams, P-lactamases have evolved at a rapid pace leading to treatment failures. Yersinia 
enterocolitica causes many gastrointestinal problems. It is an extremely heterogeneous species comprising 
more than fifty serotypes and six biotypes which differ in their ecological niches, geographical distribution 
and pathogenic potential. Though biotype 1A strains have been associated with outbreaks of Yersiniosis, 
there has been a controversy regarding their pathogenicity. The strains of Y. enterocolitica isolated from 
India belonged to biotype 1A and possessed genes for two p-lactamases namely, blaA and blaB. An earlier 
study by us reported differential expression of blaA by strains of Y. enterocolitica biotype 1A. The present 
study has been carried out to understand the molecular bases which regulate the expression of blaA in Y. 
enterocolitica biotype 1 A. We concluded that six types of blaA variants were present in strains of biotype 1 A. 
Neither amino acid substitutions in blaA nor mutations in promoter regions of WaA contributed to 
differential expression of blaA in Y. enterocolitica biotype 1 A. Rather, the secondary structures attained by 
mRNA of WaA might underlie the differential expression of blaA in Y. enterocolitica. 



Beta-lactamases and their mutated forms such as extended spectrum P-lactamases (ESBLs), AmpC (3- 
lactamases, carbapenemases and metallo P-lactamases represent a major cause of multidrug resistance 
among enteric bacteria. Following selective pressure due to widespread use of P-lactam antibiotics in 
clinical medicine, genes for P-lactamases have evolved at a very rapid pace with over 1300 unique P-lactamases 
reported in clinical isolates 1 . Yersinia enterocolitica, a member of the family Enterobacteriaceae is an important 
zoonotic pathogen. It is an extremely heterogeneous species comprising fifty serotypes and six biotypes viz. 1 A, 
IB, 2, 3, 4 and 5 which differ not only in their geographical distribution and ecological niches but also in 
pathogenic potential 2 . On the basis of pathogenicity, Y. enterocolitica can be divided into highly pathogenic 
biotype IB strains, moderately pathogenic biotype 2-5 strains and supposedly non-pathogenic biotype 1 A strains. 
Batzilla et al. 3 proposed that biotype 1 A strains represent a potential group of emerging pathogens which share 
known and putative virulence-associated features with the pathogenic bioserotypes. Studies have shown that 
strains classified as biotype 1A, in fact, represent more than one subspecies 3 . Despite the lack of classical 
determinants of virulence in biotype 1 A strains 4 these have been frequently isolated from clinical samples across 
the globe 5 " 8 . 

In India, the first food-borne outbreak due to Y. enterocolitica was reported from Tamil Nadu in 1997'. Since 
then, this pathogen has been isolated from a variety of sources - wastewater, pork, pigs, diarrheic human patients, 
ice cream, dairy products and even traditional fast food 6,9-13 . The strains of Y. enterocolitica isolated from India 611 
principally belonged to biotype 1 A. The strains of Y. enterocolitica are known to possess two P-lactamases, namely 
P-lactamase A (blaA) and P-lactamase B (blaB) 1415 . BlaA is a broad-spectrum constitutively expressed Ambler 
class A enzyme where as blaB is an inducible Ambler Class C enzyme 16 . 

The most common mode of microbial resistance to P-lactams is their enzymatic hydrolysis by P-lactamases. 
Over the time, P-lactamases have expanded their substrate spectrum with simple mutations in the gene or in the 
genetic environment of P-lactamases. Modifications in the regulatory regions, mutations in the promoter 
sequences and integration of insertion sequences containing efficient promoters have frequently been associated 
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with high-level expression of chromosomal fi-lactamases 17 . Studies 
on bla genes of Y. enterocolitica have reported differences in )3-lac- 
tamase expression and P-lactam susceptibility among different bio- 
types, focusing primarily on 'pathogenic biotypes' 1821 . The previous 
study on molecular characterization of bla genes and differences in 
PTactamase expression in biotype 1A strains was reported from our 
laboratory 14,15 . Previous studies from our laboratory also reported 
that the expression of blaB in biotype 1A strains was quite similar 
among the strains, while the expression of blah was quite variable 
among different strains 14,15 . To the best of our knowledge, till date no 
study has been carried out to understand the differential regulation of 
expression oiblaA in Y. enterocolitica biotype 1A strains. The present 
study is an extension of our previous observations and was under- 
taken to understand the molecular basis of differential expression of 
PTactamase genes of biotype 1A strains of Y. enterocolitica. This 
study also provides insights into the blah variants in biotype 1A 
strains of Y. enterocolitica. Thus, the promoter regions and the struc- 
tural regions of blah genes in biotype 1A strains of Yersinia enter- 
ocolitica isolated from India and from other parts of the world were 
PCR-amplified, sequenced and compared. Three dimensional struc- 
tures of blaA variants were modeled in-silico to co-relate the role of 
specific substitutions in differential expression of pTactamases. 

Results 

Minimal Inhibitory Concentrations of P-lactam antibiotics. Y.ente- 
rocolitica strains showed variable susceptibilities to P-lactam anti- 
biotics of different generations (Table 1). All biotype 1A strains of 
Indian and non-Indian origin except Y. enterocolitica strain E3 were 
resistant to amoxicillin and amoxicilln-clavulanate, antibiotic-inhibitor 
combination. The MIC of all biotype 1A strains to amoxicillin 
exceeded the upper range of E-test i.e > 256 mg/L. For six Indian 
strains, the MIC of amoxicilln-clavulanate combination also exceeded 
the upper range of E-test. The susceptibility of Y. enterocolitica to the 
third generation cephalosporins, cefpodoxime and cefotaxime was also 
tested by E-test. Antibiotic susceptibility of Y. enterocolitica strains to 



the third generation cephalosporin cefpodoxime varied and ranged 
from 0.125-8 mg/L. While most of the biotype 1A strains were 
either sensitive or intermediate to cefpodoxime, four strains were 
resistant to this antibiotic. Further, all Y. strains were sensitive to 
third generation cephalosporin, cefotaxime. For biotype 1A strains 
the MIC was lowest for cefotaxime and ranged between 0.032- 
0.64 mg/L. Strain E3 showed lowest MIC to all antibiotics including 
amoxicillin and amoxicillin-clavulanate. 

Specific activities of (t-lactamases. The mean specific activities of 
the blaA variants produced by Y. enterocolitica biotype 1A strains 
expressed as amount of nitrocefin hydrolyzed (umol/min/mg of 
protein) are summarized in Table 2. The specific activity of the 
variant blaA2 was quite high than the mean specific activities of 
other variants, while that of blaA4 was the least. The mean specific 
activities of the variants blaAl and blaA3 were similar to each other. 

Promoter sequences and CCDS of blaA gene. The amplification of 
blaA gene using primers blaAl 1 and blaA 12 resulted in the expected 
amplicon of 1.49 kb in 68 out of 81 strains of Y. enterocolitica. The 
band was excised from the agarose gel and sequenced after 
purification using QIAGEN purification kit. The BLAST analysis 
of the sequences confirmed these to be part of the blaA gene. 
Multiple sequence alignment of the promoter regions revealed that 
promoter sequences of all biotype 1A strains were identical (Fig. 1). 
The amplification of blaA gene using published primers 15 resulted in 
the amplification of desired product of 896 bp in all the strains. The 
band was excised from gel, purified using QIAGEN purification kit 
and sequenced. CCDSs of blaA were translated into their 
corresponding amino acid sequences and analyzed. Comparison of 
amino acid sequences of blaA using multiple sequence alignment 
revealed that though the sequence was conserved, some 
substitutions were identified. Since the initial 30 amino acids 
constitute the signal sequence and are cleaved before the mature 
enzyme is released in the periplasmic space, these were excluded 
from comparative studies. P32S and G60C was observed in 12 of 



Table 1 | MIC of selected antibiotics for V. enterocolitica biotype 1 A strains 



Strain no. 


Source 


Country 


Serotype 


AMX (mg/L) 


AMC (mg/L) 


CTX (mg/L) 


CPD (mg/L) 


V.e strain 1 


Clinical 


India 


0:6,30-6, 31 


>256(R) 


>256(R) 


0.64(S) 


1(S) 


V.e strain 2 


Clinical 


India 


0:6,30-6, 31 


>256(R) 


>256(R) 


0.125(S) 


2(S) 


V.e strain 3 


Clinical 


India 


0:6,30-6, 31 


>256(R) 


192(R) 


0.036(S) 


0.75(S) 


V.e strain 8 


Clinical 


India 


0:6,30-6, 31 


>256(R) 


128(R) 


0.125(S) 


2(S) 


V.e strain 9 


Wastewater 


India 


0:6,30-6, 31 


>256(R) 


128(R) 


0.1 9(S) 


2(S) 


V.e strain 1 0 


Wastewater 


India 


0:6,30-6, 31 


>256(R) 


192(R) 


0.38(S) 


3(1) 


V.e strain 1 1 


Wastewater 


India 


0:6,30-6, 31 


>256(R) 


64(R) 


0.75(S) 


8(R) 


V.e strain 1 2 


Clinical 


India 


0:6,30-6, 31 


>256(R) 


128 (R) 


0.1 9(S) 


8(R) 


V.e strain 19 


Clinical 


India 


0:6, 30 


>256(R) 


>256(R) 


0.125(S) 


3(1) 


V.e strain20 


Clinical 


India 


0:6, 30 


>256(R) 


>256(R) 


0.064(S) 


US] 


V.e strain 28 


Wastewater 


India 


0:15 


>256(R) 


96(R) 


0.64(S) 


US) 


V.e strain 29 


Wastewater 


India 


0:41,42 


>256(R) 


192(R) 


0.125(S) 


2(S] 


V.e strain 30 


Clinical 


India 


0:41,43 


>256(R) 


>256 (R) 


0.125(S) 


0.3(1) 


V.e strain 35 


Pig throat 


India 


ND 


>256(R) 


192(R) 


0.125(S) 


2(S] 


V.e strain 42 


Clinical 


India 


0:6, 30 


>256(R) 


192(R) 


0.25(S) 


4(1) 


y.e strain 47 


Clinical 


India 


0:6, 30 


>256(R) 


64 (R) 


0.094(S) 


2(S] 


y.e strain 48 


Clinical 


India 


0:6, 30 


>256(R) 


128(R) 


0.25(S) 


4(1) 


y.e strain 52 


Wastewater 


India 


0:6, 30-6, 31 


>256(R) 


>256(R) 


0.25(S) 


4(1) 


y.e strain 109 


Clinical 


India 


NAG 


>256(R) 


192(R) 


0.25(S) 


3(1) 


y.e strain 1 1 1 


Clinical 


India 


NAG 


>256(R) 


192(R) 


0.25(S) 


3(1) 


y.e strain E2 


Clinical 


Germany 


0:6, 30 


>256(R) 


128(R) 


0.75(S) 


8(R) 


y.e strain E3 


Clinical 


Germany 


0:6, 30 


64(S] 


0.75(S) 


0.032(S) 


0.125(S) 


y.e strain E4 


Clinical 


NK 


0:6, 30 


>256(R) 


128 (R) 


0.75(S) 


8(R) 


y.e strain E7 


Clinical 


France 


0:6, 30 


>256(R) 


1 92 (R) 


0.38 (S) 


2(S] 


y.e strain E8 


Clinical 


France 


0:6, 30 


>256(R) 


64(R) 


0.75 (S) 


4(1) 



AMX, amoxicillin; AMC, co-amoxiclav; CTX, cefotaxime CPD, cefpodoxime; 

Alphabets in parenthesis indicate antibiotic susceptibility; R-resistant, l-intermediate, S-sensitive. 

ND - not determined; NAG - non agglutinable; NK-not known. 
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Table 2 Mean specific activities of p lactamase (blaA) variants of 
Y.enterocolitica biotype 1 A strains 

blaA No. of strains in which Mean specific activity (umol/min/ 
variant the variant was present(n) mg of protein) of P-lactamase 


blaAl 12 


0.150 ± 0.02 


blaA2 1 


0.853 


blaA3 1 6 


0.153 ± 0.04 


blaA4 4 


0.030 ± 0.02 


blaA5 1 


0.10 


blaA6 1 


0.20 


All values are represented as mean ± standard 


error of mean [SEM). 



35 strains viz. Y. enterocolitica strains 12, 52, 8, 9, 42, 47, 48, 109, 1 1 1, 
E8, 10, 11; G58A in Y. enterocolitica strains E2, E4, E3, 29, 30, E7; 
R69Q in Y. enterocolitica strains E6, E2, E4, E3, 29, 30, E7; L98T and 
N101D in Y. enterocolitica strain 19. II 13V and T225S in Y. 
enterocolitica strain E4; L144V in Y. enterocolitica strains E2 and 
E4. The mutation sites in individual strains are shown in Table 3 
and summarized in Fig. 2. 

Phylogenetic analysis of blaA of Y. enterocolitica strains. The 

phylogenetic analysis of blaA of 35 strains of Y. enterocolitica 
biotype 1A was carried using the Maximum Likelihood method of 
MEGA5.5. The tree with the highest log likelihood (-939.3021) is 
shown (Fig. 3). The tree is drawn to scale, with branch lengths 



Y.enterocolitica strainl GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 201 

Y.enterocolitica strainl8 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 183 

Y.enterocolitica strainEl GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 272 

Y.enterocolitica strainl03 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 183 

Y.enterocolitica strain46 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 201 

Y.enterocolitica strainE9 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 192 

Y.enterocolitica strain7 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 231 

Y.enterocolitica strainl9 GGCAGTGGTGCTTARTCSCCATGGAGATTATAATTTCAGCCTTATG 210 

Y.enterocolitica strain20 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 208 

Y.enterocolitica strain2 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 197 

Y.enterocolitica strainl07 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 193 

Y.enterocolitica strainElO GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 196 

Y.enterocolitica strainl3 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 198 

Y.enterocolitica strainE16 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 195 

Y.enterocolitica stramE15 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 195 

Y.enterocolitica stramEll GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 199 

Y.enterocolitica strain43 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 199 

Y.enterocolitica strain3 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 200 

Y.enterocolitica strainEl2 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 192 

Y.enterocolitica strainl06 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 18 6 

Y.enterocolitica strainl7 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 191 

Y.enterocolitica strain22 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 205 

Y.enterocolitica strain4 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 198 

Y.enterocolitica strainl6 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 220 

Y.enterocolitica strain45 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 215 

Y.enterocolitica strain47 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 214 

Y.enterocolitica strain44 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 213 

Y.enterocolitica strain21 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 213 

Y.enterocolitica strainl09 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 190 

Y.enterocolitica strainllO GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 195 

Y.enterocolitica strainlll GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 195 

Y.enterocolitica strainE7 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 195 

Y.enterocolitica strainE8 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 197 

Y.enterocolitica strainl04 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 198 

Y.enterocolitica strain42 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 208 

Y.enterocolitica strain41 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 207 

Y.enterocolitica StrainE4 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 209 

Y.enterocolitica StrainE5 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 212 

Y.enterocolitica StrainE2 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 211 

Y.enterocolitica strainE3 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 210 

Y.enterocolitica Strain48 GGCAGTGGTGCTTAGTCGCCATGGAGATTATAATTTCAGCCTTATG 212 

************** ** **************************** 

Figure 1 | Multiple sequence alignment of the promoter region of MaA of representative strains of Y.enterocolitica. The transcription start site ATG, 
— 10 region and —30 regions are shown in bold face; — 10 and —30 regions were highly conserved in all the strains. 
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Table 

Sr. No 


3 Amino acid substitutions in the blaA of Y.enterocolitica biotype 1 A stra 
Strain designation Source 


ins 


Serotype 


Amino add change 


1 

1 . 


Y. enterocolitica strains 8,9, 1 0, 1 1 , 1 2,52 


^Vastewater (India) 






rJzj 


n 
z . 


Y. enterocolitica strains 42, 47,48, E8* 


Human stools (India) * Human 


(Europe) 


O-A 


r OZO 


o 
o . 


V. enterocolitica strains 109,1 1 1 


Human stools (India) 


INAL? 


roZo 


A 

4. 


Y. enterocolitica strains E2 E3 E4 E7 


Human (Europe) 




O-A Tr\ 


OOoA 


c: 
D . 


Y. enterocolitica strains 29 30 


^Vastewater (India) 






OJoA 


0. 


Y. enterocolitica strains 8, 9, 1 0, 1 1 , 1 2, 52 


Wastewater (India) 




f~VA A T 1 
U.O,oU-0,o 1 


^OU\- 


7 


Y.enterocolitica strains 42,47,48, E8* 


Human stools (India) *Human 


(burope) 


O-A TO 




Q 


Y. enterocolitica strains 109, 1 1 1 


Human stools (India) 


M AC' 




o 
y. 


Y. enterocolitica strains 29, 30 


Wastewater (India) 




\J.A \ ,AZ 




10. 


Y. enterocolitica strains E2, E3, E4, E6, E7 


Human (Europe) 




0:6,30 


R69Q 


11. 


Y. enterocolitica strain 1 9 


Wastewater (India) 




0:6,30-6,31 


L98T 


12. 


Y. enterocolitica strain 19 


Wastewater (India) 




0:6,30-6,31 


N101D 


13. 


Y. enterocolitica strain E4 


Human (Europe) 




0:6,30 


11 13V 


14. 


Y. enterocolitica strains E2, E4 


Human (Europe) 




0:6,30 


LI 44V 


15. 


Y. enterocolitica strain E4 


Human (Europe) 




0:6,30 


T225S 



represented by the number of substitutions per site. The MEGA 
analysis clustered isolates into three major groups. 

Modeling analysis. Based on the identity percentage and E-values 
(expectation values), )3-lactamase of Burkholderia multivorans (PDB 
code 3W4Q_A; Uniprot ID A9ANW2) was selected as the template 
protein due its high sequence similarity, identity and low E-value. Of 
the 20 models built using MODELLER, the 3D model with the lowest 
modeler objective function was selected as the final model. 

Protein structure validation by PROCHECK and VERIFY3D. The 

results of PROCHECK as represented in Ramachandran contour 
plot showed that ca. 90% of the residues lay in the most favoured 
regions and only 1% of residues in disallowed regions. Verify3D 
profiles which were used to identify unreliable regions that had 
been modeled improperly showed that 90% region of the protein 
model of each blaA variant scored > 0.2 which was highly significant. 

Molecular docking. Evaluation of binding affinity of pTactam 
antibiotics with class A pTactamases was performed using 
AutoDock Vina. The binding poses for each ligand molecule 
(amoxicillin, clavulanic acid) into the structures were determined 
and different poses were generated based on the total score (Dock 
score). The negative low free energy of binding of docked complexes 
and low inhibition constants indicated high affinity of pTactamase 
for these antibiotics (Table 4 & 5, Fig. 4). 

mRNA secondary structure analysis. The secondary structures of 
mRNA of the six blaA variants were analyzed. These formed 
branched secondary structures. The structure of blaA2 showed 
very high level of branching, with a number of bulge loops, interior 
loops and a smaller region of double stranded RNA helix. On the 
other hand the mRNA secondary structure of blaA4 was mostly 
composed of a double-stranded RNA helix of stacked base pairs 
with little branching (Fig. 5). 



Discussion 

The MIC of different PTactam antibiotics for Y. enterocolitica bio- 
type 1A strains using E-test revealed that all strains were resistant to 
amoxicillin and amoxicillin-clavulanate. Y. enterocolitica strains 
have been reported to be 'naturally resistant' to ampicillin 30,31 . 
Although the MIC of amoxicillin-clavulanate was perceptibly less 
than that of amoxicillin, it did not change the resistance phenotype, 
implying the existence of inhibitor resistant blaAs in Y. enterocolitica 
biotype 1A. Earlier studies have shown that Y. enterocolitica strains 
of serotypes O: 3, O: 5 and O: 9 31 were all resistant to ampicillin and a 
few were resistant to amoxicillin-clavulanate. Szych et al. 32 reported 
that though all strains of bioserotypes 4/0: 3 and lB/O: 8 were 
resistant to ampicillin, yet 90% of these were sensitive to amoxicil- 
lin-clavulanate. They also reported that strains of bioserotype 2/0: 9 
were resistant to both ampicillin and amoxicillin-clavulanate. Recent 
studies 8,33 also showed that strains of Y.enterocolitica bioserotype 4/ 
O: 3 were sensitive to amoxicillin-clavulanate where as that of bio- 
serotypes 2/0: 5, 27 and 2/0: 9 were resistant. From the available 
information it is apparent that both the biotype and serotype play 
role in affecting the susceptibility of Y. enterocolitica to amoxicillin- 
clavulanate. Since all strains in the present study belonged to biotype 
1A, it is not possible to comment on the influence of biotype on 
antibiotic susceptibility; but based on our findings it may be con- 
cluded that serotype had no effect on the susceptibility of Y. enter- 
ocoloitica biotype 1 A strains to amoxicillin-clavulanate. The biotype 
1 A strains were either sensitive or showed intermediate susceptibility 
to the third generation cephalosporin, cefpodoxime and only a few 
strains were resistant. This is the same as observed in an earlier 
study 8 . Y. enterocolitica strain E3 was sensitive to all the PTactam 
antibiotics tested and showed lowest MIC to all the PTactam anti- 
biotics tested including amoxicillin and amoxicillin-clavulanate 
combination. 

The mean specific activities of variants blaAl and blaA3 were 
almost similar and the mutations/substitutions which lead to these 
variants had no effect on the enzyme activity. Similarly mutations 
that generated variants blaA5 and blaA6 did not affect the enzyme 



. .31 


. . 58 


. 60... 


. 69 


. 98... 


. . 101 


113 


144 


225 


aminoacid position 


. . P. . . . 


. . . G... 


. G. . 


.R. . . 


. .K. . 


. . . N . . . . 


. . . . I 


L 


T 


blaA3 (n=16) 


. .S. . . . 


. . . G... 


.C. . 


. R. . . 


. . K. . 


. . . N . . . . 


. . . . I 


L 


T 


blaAl (n=12) 


. .P. . . . 


. . .A... 


.G. . 


■ Q- ■ ■ 


. . K. . 


. . . N . . . . 


. . . . I 


L 


T 


blaA4 (n=4) 


. .P. . . . 


. . . G... 


.G. . 


. R. . . 


. . T . . 


. . . D. . . . 


. . . . I 


L 


T 


blaA2 (n=l) 


. .P. . . . 


. . .A... 


. G. . 


■ Q- ■ ■ 


. . K. . 


. . . N . . . . 


. . . . V 


V 


S 


blaA5 (n=l) 


. .P. . . . 


. . .A... 


. G. . 


.Q. . . 


. . K. . 


. . . N . . . . 


. . . . I 


V 


T 


blaA6 (n=l) 



Figure 2 | The blaA variants present in biotype 1 A strains of Y.enterocolitica, where n represents the number of strains in which this type of amino acid 
sequence was observed. 
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Group I 



■ Y. e strain 1 9 



Y. e strain 106 
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Figure 3 | Phylogenetic analysis of MaA gene of Y.enterocolitica biotype 
1A strains. 



activity significantly. However, the mutations leading to variants 
blaA2 and blaA4 significantly increased or decreased the specific 
activities of the enzyme respectively. 

Bacteria have been reported to hyper produce |3-lactamases by 
different mechanisms viz. mutation in the promoters, integration 
of insertion sequences containing efficient promoters or point 
mutations in the structural gene 17 . In the family Enterobacteria- 
ceae, mutations and insertions, mostly in the promoter regions of 
the (3-lactamase genes, have been associated with hyper production 
of P-lactamases 34 36 . In the present study, we attempted to under- 
stand the molecular basis of differential expression of chromosomal 
blaA of biotypel A strains of Y. enterocolitica. For this, the nucleotide 
sequences of promoter regions of blaA of biotype 1 A strains were 
PCR-amplified and analyzed. Genetic events that affect the tran- 
scription efficacy of bacterial promoters usually occur in the — 35 
or — 10 regions or sometimes in the space between these regions. It is 
well known that promoter strength is linked to the sequence of the 
— 35 and —10 boxes, creating the 70 bp contact region for RNA 
polymerase. The —10 sequence TATAAT and the —35 sequence 
GTGGTG of blaA of Y. enteocolitica were similar to the canonical 
E.coli sequence indicating these were strong promoters. Alignment 
and comparison of blaA of biotype 1 A strains by Clustal revealed 
that the ribosomal binding site (RBS), the —10 region and —35 
regions were identical in all strains indicating that the promoter 
region had no role in differential expression of P-lactamases. Thus 
it was argued that the differential expression of class A beta-lacta- 
mase enzyme in Y. enterocolitica biotype 1A might be due to the 
changes in gene sequence of blaA. 

To investigate the role of changes in gene sequence of blaA in 
differential expression of enzyme, blaA genes of biotype 1 A strains 
were PCR-amplified and their CCDSs translated into corresponding 
amino acid sequences. The four significant motifs conserved across 
the Ambler class A P-lactamase viz. SXXK, SDN, EXXLN and KTG, 
(where X represents any amino acid) were present in blaA of all Y. 
enterocolitica strains, their respective sequence and position being 

75 STpK 78 ; 136 SDN 13 8] 171 EPDLN 175 and Qf mg four 

conserved domains 171 EPDLN 175 was located in the omega loop, a 
structural domain constituting a part of active-site pocket. Clustal 
alignment and comparison of amino acid sequences of structural P- 
lactamase blaA of biotype 1 A strains revealed that amino acid sub- 
stitutions were present at some places. Since the initial 30 amino 
acids constitute the signal sequence and are cleaved before the 



Table 4 


Estimated inhibition constants, free energy c 


f binding, H-bond interactions and hydrophobic interactions between amoxicillin and 


blaA variants of Y. enterocolitica biotype 1 A 






blaA 


Estimated Inhibition 


Free Energy of 






variant 


Constant, Ki 


Binding (kcal/mol' 




Interacting Residues 


blaAl 


57. 1 2 uM 


-5.79 


Hbonds 
Hydrophobic 
interactions 


75S, 135S, 137N, 175N, 242S 

78K, 109S, 1 10Y, 171 E, 172P, 22 IT, 243G, 244D, 279R 


blaA2 


14.30 uM 


-6.61 


Hbonds 
Hydrophobic 
interactions 


75S, 78K, 1 35S, 22 IT, 242S, 279R 

74C, 109S, HOY, 137N, 171 E, 172P, 175N, 240T, 241 G 


blaA3 


4.31 uM 


-7.32 


Hbonds 
Hydrophobic 
interactions 


75S, 109S, 135S, 137N, 171,242,279 
78K, 1 1 0Y, 1 34Y, 1 75N, 22 1 T, 240T 


blaA4 


21.79 uM 


-6.36 


Hbonds 
Hydrophobic 
interactions 


1 35S, 1 37N, 1 75N, 240T, 242S 

75S, 78K, 109S, 1 10Y, 171E, 22 IT, 241 G, 243G 


blaA5 


1 16.57 uM 


-5.37 


Hbonds 
Hydrophobic 
interactions 


75S, 109S, 135S, 137N, 175N, 242S 
78K, 1 10Y, 171 E, 172P, 22 IT, 243G 


blaA6 


15.17 uM 


-6.57 


Hbonds 
Hydrophobic 
interactions 


75S, 1 35S, 1 37N, 1 71 H, 242S, 279H 
78K, HOY, 175N, 22 IT, 240T 
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Table 5 


Estimated inhibition constants, free energy of binding, H-bond interactions and hydrophobic interactions between clavulanic acid 


and blaA variants of Y. enterocolitica biotype 1 A 






blaA 


Estimated Inhibition 


Free Energy of Binding 






variant 


Constant, Ki 


(kcal/mol) 




Interacting Residues 


II A 1 

blaAl 


812.13 uM 


-4.22 


Hbonds 
Hydrophobic 
interactions 


75S, 1 35S, 242S 

78K, HOY, 137N, 22 IT, 239K, 240T, 24 1G, 249N 


LI AO 

blaA/! 


it in k A 
OAV.VO uM 


A O C 

-4.35 


Hbonds 
Hydrophobic 
interactions 


75S, 78K, 1 35S 

HOY, 137N, 175N, 240T, 24 1G, 242S 


LI AO 

blaAJ 


o on "7n k A 
JVZ./2 uM 


-4.65 


Hbonds 
Hydrophobic 
interactions 


75S, 135S, 22 IT, 239K, 242S 

78K, HOY, 134Y, 137N, 240T, 24 1G 


blaA4 


1.11 uM 


-4.03 


Hbonds 
Hydrophobic 
interactions 


75S, 1 35S, 1 37N 

74C, 1 09S, 1 1 0Y, 1 71 E, 1 75N, 242S 


blaA5 


159.42 uM 


-4.98 


Hbonds 
Hydrophobic 
interactions 


135S, 240T, 242S, 239K 

75S, 221T, 241G, 248T, 249N, 277P, 279R, 282L 


blaA6 


150.09 uM 


-5.18 


Hbonds 
Hydrophobic 
interactions 


1 35S, 240T, 242S 

75S, 22 IT, 239K, 241 G, 248T, 249N, 277P, 279R, 282L 



mature enzyme is released in the periplasmic space, these were 
excluded from the comparative studies. The most frequent point 
substitutions observed in 12 of the 35 Y. enterocolitica strains were 
a substitution from P32S and G60C. As proline and serine are both 
polar amino acids, substitution of one with the other will affect only 
the bulkiness of the side chain of the protein. Glycine is a simple 
nonpolar amino acid its replacement with cysteine, a polar sulfur 
containing amino acid may affect the properties of the side chain of 



the protein. The substitution of glycine with alanine G58A in Y. 
enterocolitica strains would also affect only the bulkiness of the side 
chain because both are nonpolar aliphatic amino acids. As arginine is 
a positively charged amino acid, while glutamine is a polar uncharged 
amino acid, substitution of R69N in Y. enterocolitica strains may also 
affect enzyme activity. Lysine is a positively charged amino acid while 
threonine is a polar uncharged amino acid; substitution K98T is 
expected to affect the side chain of the protein. Asparagine is a polar 




Figure 4 | Molecular docking analysis of the six blaA variants with amoxicillin and clavulanic acid. Molecular interactions of docked blaAl (a & b), 
blaA2 (c & d), blaA3 (e & f), blaA4 (g & h), blaA5 (i & j), blaA5 (k & 1) and blaA6 (m & n) with amoxicillin and clavulanic acid respectively. The two 
antibiotics are represented by stick and the interacting amino acids are shown in dark color. 
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Figure 5 | The mRNA secondary structures of blaAl (a), blaA2 (b), blaA3(c), blaA4 (d), blaA5 (e) and blaA6 (f) predicted by mfold based on complete 
coding sequences of different blaA variants. 



uncharged amino acid, while aspartate is a negatively charged amino 
acid; substitution of N101D in Y. enterocolitica strain 19 may influ- 
ence enzyme activity due to change in the properties of the side chain. 
Substitution of I111V and T225S in Y. enterocolitica strain E4 and 
L144V in Y. enterocolitica strains E2 and E4 would also affect only the 
bulkiness of the side chain of the protein without any profound effect 
on enzyme activity because of the similarity in the chemical nature of 
the substituent and the substitute. For example isoleucine and valine 
are nonaliphatic amino acids, threonine and serine are polar amino 
acids and leucine and valine both nonpolar amino acids. 

Phylogenetic analysis based on amino acid sequences of blaA clus- 
tered the biotype 1 A strains into 3 major groups. In our earlier study 
also the Y. enterocolitica biotype 1A strains clustered into three 
groups based on the profiles generated by the restriction analysis 
ofWaA 15 . 

The analysis of amino acid sequences of blaA of Y. enterocolitica 
biotype 1A strains showed that although none of the substitution/ 
mutation occurred in the four major motifs conserved across class A 
P-lactamases, or in the omega-loop region of blaA, yet there were 
differences in the specific activities of blaA. Amino acid substitutions 
which do not occur in the active site of the enzyme may create diverse 
local changes in the 3D structure of the enzyme and increase its 
conformational flexibility which may in turn affect the binding affin- 
ity with the substrate. To evaluate the role of such amino acid sub- 
stitutions on binding affinity we classified the blaA of biotype 1A 
strains of Y. enterocolitica into six variants based on the mutation/ 
substitution in the amino acid sequence. Three dimensional model- 
ing of the blaA variants was carried out to co-relate structural 
changes with the higher enzyme specific activities and MIC of dif- 
ferent strains. In-silico docking of the blaA variants was carried out 
with amoxicillin and clavulanic acid to analyze if amino acid sub- 
stitutions affected the residues involved in hydrogen bonding and 
hydrophobic interactions in the enzyme-ligand complex. It is well 
known that hydrogen bonds and hydrophobic interactions not only 
play a crucial role in molecular recognition and the overall stability of 
the protein structure; these also stabilize the protein-ligand com- 
plexes. Amoxicillin and clavulanic acid were chosen for docking 
experiments as both are penicillins and blaA is a penicillinase. 
Intermolecular hydrogen bonds formed between the modeled struc- 
tures of blaA variants and amoxicillin/clavulanic acid were deter- 
mined. Pymol analysis revealed that each blaA variant had its own 
drug-interacting active-site amino acids. In all blaA variants S-135, 
A- 137 and S-242 participated in hydrogen bonding with amoxicillin, 
and S-135 in hydrogen bonding with clavulanic acid. Kumar et al. 37 
reported that the number of the residues involved in hydrogen bond 



interactions was an indicator of binding affinity of the enzyme with 
the ligand which might be decreased in mutant forms of class A (3- 
lactamase. However, in our study no co-relation was observed 
between the number of residues involved in hydrogen bond interac- 
tions in blaA variants and their MIC for amoxicillin or clavulanic 
acid. Yi et al. 38 also observed that the cetadazidime MIC and the 
relative frequency of substitutions in PenA, a class A p-lactamase 
of Burkholderia thailandensis did not co-relate well. It may be worth- 
while to mention that the Y. enterocolitica isolates used in the present 
work were isolated from environmental as well as clinical settings. 
Variations/mutations associated with their respective blaAs might 
have evolved over time either in nature or in clinical settings in 
patients treated with antibiotics. Salverda et al. 39 have also reported 
that most of the mutations that evolved in P-lactamases in clinical 
settings do not clearly co-relate with the stability or enhanced hydro- 
lytic capability of the enzyme. 

Variation in amino acid sequence is not the only factor that might 
contribute to differences in specific activities of P-lactamases. 
Enzyme activity may also be affected at the level of translation. In 
this regard, the secondary structure of mRNA plays an important 
role. When the secondary structures of mRNA of blaA variants were 
predicted using the web server mfold it was found that except for 
blaA4, all variants formed branched secondary structures. The 
mRNA secondary structure of blaA2 showed extensive branching 
and a lesser region of double stranded RNA helix compared to other 
blaA variants. On the other hand the mRNA secondary structure of 
blaA4 comprised mainly of double-stranded RNA helix of stacked 
base pairs with little branching. Several studies have shown that 
longer and stronger-paired zones in mRNA secondary structure tend 
to be more stable, preventing the ribosome from breaking the pairing 
and helping translation to proceed 40,41 . A lower translation rate 
means lesser number of copies of protein and hence a lower enzyme 
specific activity. Thus the secondary structures attained by mRNA of 
blaA2 and blaA4 might be the probable reasons for the higher and 
lower specific activities of their respective P-lactamases. 

Overall the present work revealed that variations in promoter 
regions were not responsible for differential expression of blaA in 
biotype 1A strains of Y. enterocolitica. Analysis of amino acid 
sequences revealed that six types of blaA variants were present among 
biotype 1A strains of Y. enterocolitica. In-silico analysis of the blaA 
variants showed that amino acid substitutions did not contribute to 
the differential expression of the enzyme activities either. Rather the 
secondary structures attained by mRNA of blaA might underlie the 
differential expression of blaA in Y. enterocolitica. Further studies on 
secondary structures of mRNA of blaA gene will be helpful in better 
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understanding the mechanisms that might underlie the differential 
expression of P-lactamases in Y. enterocolitica. 

Methods 

Bacterial strains. A total of 81 strains of Y. enterocolitica belonging to biotype 1A 
maintained in our laboratory (University of Delhi South Campus) were examined. Of 
these, sixty-five strains were isolated from India, which comprised 35 clinical, 18 
wastewater, 7 swine and 5 pork isolates. All isolates were authenticated by, and have 
been deposited with Yersinia National Reference Laboratory and WHO 
Collaborating Center, Pasteur Institute, Paris (France). All these strains have also 
been deposited with the National repository i.e. Microbial Type Culture Collection 
(MTCC) and Gene Bank located at Institute of Microbial Technology, Chandigarh. 
Ten strains of Y. enterocolitica biotype 1A belonging to European and American 
origin were obtained from E. Carniel (Yersinia National Reference Laboratory and 
WHO Collaborating Center, Pasteur Institute, Paris, France). Another six biotype 1A 
strains of European origin were provided by J. Heesemann (Max von Pattenkofer 
Institute, Munich, Germany). The serovars, source of isolation, country of origin and 
reference laboratory accession numbers of these strains have been reported 
previously 22 . All strains were maintained on trypticase soy agar at 4°C. 

Determination of the Minimal Inhibitory Concentration (MIC). MIC of 

amoxicillin, amoxicillin -clavulanate, cefotaxime and cefpodoxime was determined 
for 25 strains of Y. enterocolitica biotype 1 A using E-test (bioMerieux Inc., MO, USA) 
following the manufacturer's instructions. Briefly, each strain was grown on Muller- 
Hinton agar plates at 25°C for 18-24 h. Single colonies from the agar plates were 
suspended in 1 ml of normal saline until a turbidity of 0.5 McFarland standard was 
reached (equivalent to ca.1-2 X 10 8 CFU/ml). Sterile glass spreader was used to 
spread 100 ul of the cell suspension on the surface of Muller-Hinton agar plates. The 
E-test strips were placed on the surface of inoculated Muller-Hinton agar plates and 
incubated at 28°C for 16-18 h. The MIC values were recorded as the minimum 
concentration of the antibiotic (ug/ml) where no visible growth of the test organism 
was observed and were interpreted according to the guidelines of Clinical and 
Laboratory Standards Institute 23 . 

Spectrophotometeric assay of p-lactamases. Thirty-five strains of Y. enterocolitica 
biotype 1A were used for assaying |3-lactamase activity as reported by Stock and 
Weidmann 18 with slight modification. A loopful of bacterial growth was inoculated in 
10 ml of LB broth and incubated overnight at 28°C. The bacteria were sedimented by 
centrifugation at 4 U C, washed thrice with chilled phosphate buffer (0.05 M, pH 7.0) 
and finally resuspended (1 g wet weight) in 5 ml of sonication buffer. The cells were 
disrupted by sonication on ice with three pulses of 30 sec each. The clear cell lysate, 
prepared by centrifugation at 12,000 g for 30 min at 4°C was used to assay the |3- 
lactamase activity. The protein content of the lysates was estimated by Bradford's 
method 24 , ^-lactamase activity was assessed spectrophotometrically by hydrolysis of 
nitrocefin. The assay mixture contained 1 10 ul nitrocefm (0.5 ug/ul), 30 ug of crude 
cell lysates containing ^-lactamase in a final volume of 1.0 ml of 50 inM phosphate 
buffer (pH 7.0). [3-lactamase activity was monitored by measuring decrease in the 
optical density at 486 nm for 10 min at 28°C. The specific activity of the enzyme was 
expressed as umol of nitrocefin hydrolyzed/min/mg of protein and was calculated 
based on the molar extinction coefficient of 20,500 M _1 cm _1 for nitrocefin 25 . 

Preparation of genomic DNA. Each strain was grown overnight in trypticase soy 
broth at 28°C. One ml of the culture was pelleted at 8,000 X g for 10 min. The 
genomic DNA was prepared from each strain by using DNeasy Tissue kit (Qiagen, 
Hilden, Germany) with modifications for Gram-negative bacteria according to the 
manufacturer's recommendations. Purified DNA was eluted in sterile water and 
quantitated spectrophotometrically at 260 nm. 

Primer designing for PCR amplification of blaA gene including the promoter and 
downstream region. For amplification of partial coding sequence of blaA including 
its promoter, primers Allf 5 ' CAAATCGCCAGCAGAACAGA 3' and Al2r 
5 ' GGGCAGCGAGCAGGGGTAAA 3' were designed by retrieving the whole 
genome sequence of Y. enterocolitica biotype IB and the corresponding sequence of 
blaA along with the 350 bp each upstream and downstream of the gene. The primers 
were synthesized from Microsynth GmbH (Balgach, Germany). 

PCR amplification of the promoter region and partial coding sequence of blaA was 
carried out for 81 strains of Y. enterocolitica in a Thermal Cycler (My Cycler™, Bio- 
Rad, CA, USA). The PCR reaction mixture comprised of 1 X PCR buffer ( 10 mM Tris- 
HCL, 1.5 mM MgCl 2 , 1.5 mM KC1 and 0.1% Triton X-100), 200 uM of each of the 
four dNTPs (MBI Fermentas GmbH, Germany), 10 pmol each of forward and reverse 
primers, 2 U of Taq DNA polymerse (NewEngland Biolabs, MA, USA) and 50-100 ng 
of genomic DNA in a total volume of 25 ul. For PCR, the initial denaturation was 
performed at 95 C for 5 min followed by 25 consecutive cycles of 30 sec at 95 C, 
30 sec at 56"C and 90 sec at 72°C and a final extension for 10 min at 72°C. The PCR 
products were analyzed by electrophoresis on 1% (w/v) agarose gel. The gels were 
stained with ethidium bromide (0.5 |J.g/ml) and visualized under UV transilluminator. 

PCR amplification of blaA gene. PCR amplification of complete coding sequence 
(CCDS) of blaA was performed for 35 strains using published primers 15 in a Thermal 
Cycler (My Cycler™, Bio-Rad, CA, USA). The PCR conditions were same as that for 
partial coding region of blaA and its flanking regions. 



Sequencing of CCDS of blaA and its upstream and downstream regions. PCR 

amplicons were excised from gel, purified using QIA Quick Gel Extraction Kit (Qiagen, 
Hilden, Germany). The purified products were sequenced using the Big Dye 
Terminator Cycle Sequencing Ready Reaction kit and analyzed by capillary 
electrophoresis on an ABI PRISM 310 Genetic Analyzer (Applied Biosystems, CA, 
U.S.A). The nucleotide sequences amplified using published primers 15 were aligned 
with those sequences which were amplified using primer pair Allf and A12r. The 
sequences which were common to the amplicons obtained using the two sets of primers 
and the extra sequences which were not common were used to construct the full 
sequence of blaA including the upstream and downstream sequences. CCDSs of blaA 
of 35 strains of Y. enterocolitica biotype 1A were translated into their corresponding 
amino acid sequences using the software expasy (www.expasy.org/translate). The 
amino acid sequences of blaA were aligned by Clustal . (www.ebi.uc.ak/clustal ). 

Phylogenetic analysis of blaA of Y. enterocolitica strains. The amino acid sequences 
of blaA of 35 strains of Y. enterocolitica strains biotype 1A were also analyzed from 
evolutionary perspectives using the software MEGA5.2 (Molecular Evolutionary 
Genetic Analysis). The evolutionary history was inferred by using the Maximum 
Likelihood method based on the JTT matrix-based model (www.megasoftware.net). 

Selection of Y. enterocolitica strains for modeling of blaA. To investigate the role of 
point substitutions in amino acid sequence in differential expression of blaA, analysis 
of protein models of blaA variants was carried out. After analysis by Clustal, amino 
acid sequences of blaA of Y. enterocolitica strains representing type/(s) of 
substitution/ (s) were selected for modeling of their P-lactamases. The blaA of Y. 
strains 8, 19, 20, 29, E2 and E4 were selected for modeling and accordingly named as 
blaAl, blaA2, blaA3, blaA4, blaA5 and blaA6. 

Homology modeling of p-lactamases. Since the three dimensional (3D) structure of 
blaA of Y. enterocolitica is not known, the 3D structures of ^-lactamase variants 
blaAl, blaA2, blaA3, blaA4, blaA5 and blaA6 were predicted by homology modeling. 
The pairwise alignment between the target and template sequences was performed 
with PDB-BLAST. Three-dimensional structures of blaAl, blaA2, blaA3, blaA4, 
blaA5 and blaA6 were built using the software MODELLER 9.12 (http://salilab.org/ 
modeller/). Twenty models satisfying the spatial restraints based on target- template 
alignment were built for each protein, of which the 3D model with the lowest modeler 
objective function was selected. Following the modeling process, the structure was 
validated by PROCHECK and VERIFY3D. The PROCHECK program performs a 
detailed analysis of the stereochemistry of the modeled protein structure in the form 
of Ramachandran plot by plotting the phi-psi torsion angles. Population of phi-psi 
torsion angles in the plot gives an assessment of the overall quality of the modeled 
structure and also highlights regions that may need further investigation 26 . The 
quality of the three-dimensional models was also assessed by VERIFY3D, which tests 
the compatibility of a protein structure with its amino acid sequence. For high- 
resolution, experimentally determined structures, the VERIFY3D scores are positive 
and consistently high i.e > 0.2 27 . 

Molecular docking. Molecular docking analysis of the modeled protein structures was 
carried out to identify if substitutions in the amino acid sequences of blaA resulted in 
variations in their binding affinity to [3-lactam antibiotic amoxicillin and the [3- 
lactamase inhibitor clavulanic acid. The structures of amoxicillin and clavulanic acid 
were retrieved from the PubChem compound database (http://pubchem.ncbi.nlm.nih. 
gov/) 28 . Subsequently the downloaded SDF file was converted to PDB format using 
iBabel. Protein structures were pre-processed for docking using latest AutoDockTools 
(ADT) version by adding polar hydrogen atoms, removing all non-protein molecules 
from the modeled structure, adding Kollman charges to the protein and converting it to 
PDBQT format. Ligand structures were also prepared using ADT, using the default 
method for preparing ligands for docking that added hydrogen and charges. The 
default rotatable bonds were accepted as well, and the structure was converted to 
PDBQT format The search space for docking was determined visually by centering the 
Grid Box to the known binding site of the ligand and expanding the dimensions of the 
cubic search space to completely encompass the site. Docking was performed using 
AutoDock Vina (http://vina.scripps.edu/index.html) with default parameter settings 
other than the search space specification described above, and the mean predicted 
binding affinity from the set of predicted binding poses was accepted as the true binding 
affinity for each docking run. The maximum number of poses per ligand was set to 10. 
The predicted binding of the ligand was compared to the known binding sites of |3- 
lactamase in order to make predictions about differential responses with respect to each 
of the drugs. Interaction between ligand and target protein was analyzed by PyMOL 29 . 

Analysis of mRN A secondary structure. The secondary structure of mRNA of blaA 
variants was predicted using the mfold webserver at default parameters (http://mfold. 
rna.albany.edu/). The mfold predicts energetically optimal secondary structure of an 
RNA molecule taking into account many realistic physical parameters that affect the 
RNA folding, such as pH, temperature and local biases of RNA. 
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